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I .  INTRODUCTION 


A.  background 

The  Ballistic  Research  Laboratory  is  currently  in  the  process  of  acquir¬ 
ing  a  target  enclosure  to  facilitate  destructive  terminal  ballistic  testing  of 
chemical  explosives  (CE) ,  armor  and  kinetic  energy  (KE)  penetrators  by  safe 
containment  of  blast,  fragments  and  resultant  combustion  products.  The 
present  investigation  is  based  on  a  preliminary  concept  of  the  firing  range 
as  shown  in  Figure  1.  The  target  is  located  inside  the  hemispherical 
enclosure  at  the  end  of  a  long  concrete  pipe-guide.  The  gun-launched  projec¬ 
tile  travels  through  the  pipe-guide  and  enters  the  enclosure  through  a  .914  m 
diameter  hole.  The  target  interaction  with  the  projectile  is  monitored 
photographically  with  flash  X-ray  equipment,  and  penetration  velocity  is  obtain¬ 
ed  using  velocity  screens  and  electronic  counters.  An  air  exhaust  system 
mounted  at  the  rear  of  the  structure  operates  during  the  test  and  draws  back 
aerosolized  material  out  of  the  enclosure  after  a  test  and  traps  it  in  filters 
in  the  exhaust  ducting.  A  large  sliding  door  with  a  configuration  to  match 
the  curvature  of  the  hemispherical  wall  allows  equipment  access  inside  the 
enclosure.  The  door  is  sealed  to  the  wall  with  a  pressurized  hose  seal  along 
its  perimeter.  The  entire  structure  is  built  to  contain  blast  and  fragments, 
to  trap  aerosolized  materials  and  to  permit  photographic  observation  of  the 
test . 


A  significant  problem  associated  with  the  enclosed  range  tests  is  the 
overpressure  resulting  from  shock  loading  as  well  as  rapid  heating  of  the  air 
within  the  enclosure  as  the  penetrator  and  the  target  are  torn  apart  during 
their  encounter  as  shown  by  R.  Abrahams  et  al.1  The  structure  must  survive 
both  the  reflected  and  the  residual  overpressures  induced  by  the  interaction 
until  ambient  conditions  are  reached  due  to  venting  out  to  the  atmosphere 
through  the  exhaust  system. 

B.  Objectives 

Since  the  key  element  of  the  AHKELS  (Advanced  High  Kinetic  Energy  Launch 
System)  range  is  the  enclosure  structure.  The  Blast  Dynamics  Branch  was 
assigned  to  estimate  the  overpressure  loading  on  the  wall  and  analyze 
dynamic  response  of  the  preliminary  configuration  at  critical  locations  and 
assure  structural  integrity  from  a  conservative  viewpoint.  The  choice  of  a 
hemispherical  configuration  was  influenced  by  an  earlier  investigation  by 
N.  J.  Huffington  et  al,2  who  demonstrated  the  effectiveness  of  such  a 
protective  structure. 

In  the  absence  of  any  available  experimental  data, it  was  decided  to 
obtain  a  theoretical  estimate  of  the  transient  and  residual  overpressure 
loading  due  to  a  centrally  located  equivalent  charge  weight  at  the  base.  The 


1r.  Abrahams,  R.  Peterson,  and  B.  Bertrand,  " 'Measurement  of  Blast  Pressure 
Produced  by  Impact  of  Kinetic  Energy  Penetrator  on  a  Steel  Target,"  BRL 
Memorandum  Report  ARBRL-MR-0298Z,  January  1980  (ADB  045141L). 

2 

N.  J.  Huffington  and  S.  R.  Robertson,  " Containment  Structures  Versus 
Suppressive  Structures,"  bRL  Metro  randum  Report  BRL-MR-  25  i  7,  February  1976 
(ADA  021973 }.  _ 


subsequent  objective  was  to  perform  an  approximate  conservative  static 
analysis  for  an  initial  estimate  of  wall  thickness.  Finally  the  dynamic, 
elasto-plastic,  large  deflection  response  of  the  shell  configuration  clamped 
to  a  horizontal  rigid  foundation  was  studied  to  indicate  critical  locations 
where  peak  strains  or  deflections  could  occur. 


II.  ESTIMATION  OF  TRANSIENT  LOADS 


The  transient  loads  were  estimated  under  the  assumption  that  the  test 
firing  of  penetrator  rounds  would  generate  overpressures  inside  the  contain¬ 
ment  chamber  similar  to  those  caused  by  an  internal  blast  due  to  an  equivalent 
central  charge  weight  of  29.03  Kg*  at  the  base  as  shown  in  Figure  2.  Assuming 
the  walls  to  be  rigid,  the  symmetry  of  the  charge  and  the  structure  about  a 
vertical  axis  through  the  center  indicates  uniform  distribution  of  internal 
reflected  loading  upon  the  structure.  For  the  estimation  of  peak  reflected 
overpressure,  a  conservative  cube-root  scaling  law3  is  employed  to  compute  the 
scaled  distance,  Z,  of  the  wall  from  the  charge  location  in  the  form 


Z  = 


(1) 


where  is  the  equivalent  charge  weight  and  R  is  the  distance  of  the  wall 
from  the  charge  location . 

Once  the  scaled  distance  is  known  the  reflected  parameters  such  as  peak 
overpressure,  impulse,  time  of  arrival  and  duration  time  of  the  shock  loading 
could  be  estimated  from  compiled  air  blast  tables. 4, 5  The  decay  of  the  reflec 
ted  overpressure  is  assumed  to  obey  the  modified  Friedlander  exponential  decay 
equation  which  can  be  written  as 


-at/t 

P  =  P  [1  -  t/t  ]e  °  (2) 

r  mL  oJ 

where  tQ  is  the  positive  phase  duration  of  the  impulse,  Pm  is  the  peak  reflec¬ 
ted  overpressure  and  t  is  the  elapsed  time  from  arrival  of  blast  wave  at  inner 
surface  of  the  hemisphere.  The  exponential  decay  parameter,  ot,  is  calculated 
from 


As  per  equivalent  charge  data  provided  by  L.  Giglio-Tos ,  Armor  Mechanics 
Branch ,  Terminal  Ballistics  Division,  BRL,  September  1980. 

3 

Engineering  Design  Handbook ,  ", Explosions  in  Air,  Part  One,"  AMC  Pamphlet 
AMCP  706-181,  Headquarters ,  US  Army  Materiel  Command ,  pp.  3-5,  July  1974. 

^ B .  Soroka,  "Air  Blast  Tables  for  Spherical  50/50  Pentolite  Charges  at  Side-On 
and  Hormal  Incidence,"  BRL  Memorandum  Report  ARBRL-MR-02975y  December  1979 
(ADA  080537). 

H .  Goodman,  "Compiled  Free  Air  Blast  Data  on  Bare  Spherical  Pentolite,"  BRL 
Report  BRL-R-1092,  February  1960  (ADA  235278) • 
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Figure  2.  Sectional  View  Through  the  Hemispherical  Containment  Structure. 


where  I  is  the  reflected  impulse.  Values  of  P  ,  I  ,  and  t  were  taken  from 
r  r  r  r  o 

tables  listed  in  Reference  4.  When  a  is  determine^  the  complete  reflected 

pressure  versus  time  loading  on  the  structure  can  be  predicted. 


III.  ESTIMATION  OF  QUASI-STATIC  LOADS 

Quasi-static  pressures  immediately  following  the  reflected  pressure 
were  predicted  assuming  that  the  heat  of  combustion  of  TNT  is  used  totally  to 
heat  the  air  within  the  enclosure6  A  relationship  for  the  resultant  increase 
in  pressure  is 

0.4  hcW 

AP  =  - y— ^  ,  kPa,  (4) 

where 

V  =  1513.9  m3,  the  internal  volume  of  the  enclosure, 

Wg  =  29.03  Kg,  weight  of  the  explosive  charge,  and 

hc  =  13.5  kJ/g,  the  heat  of  combustion  of  TNT. 

An  internally  pressurized  structure  vents  the  pressure  to  the  surroundings 
through  openings  in  its  walls,  causing  a  slow  decay  to  ambient  conditions 
as  shown  by  Kinney  and  Sewell7  and  is  computed  from 

Ln  P  =  Ln  P  -  .315  (A  /V)  t  ,  (5) 

o  v  s 

where 

t  =  elapsed  time  in  ms 

P  =  absolute  pressure  at  t 

A  =  2.33  m2,  the  available  vent  area, 
v 


6 

Eduard  M.  Weyer ,  Editor-in-chief,  Annals  of  the  New  York  Academy  of  Sciences, 
Vol.  152,  " Prevention  of  and  Protection  Against  Accidental  Explosion  of 
Munitions,  Fuels  and  other  Hazardous  Mixtures Published  by  the  Academy, 

2  East  Sixty-Third  Street,  New  York,  NY  10021 ,  p.  317. 

7 

G.  F.  Kinney  and  R .  G.  S.  Sewell,  "Venting  of  Explosives,"  NWC  Technical 
Memorandum  No.  2448,  July  1974. 


The  long-term  duration  of  the  decay  is  essentially  due  to  the  relatively  small 
vent  area  available, causing  a  slow  pressure  decay  to  the  atmosphere. 

The  blow-down  time,  t  ,  required  to  reduce  the  residual  overpressure  to 

8  . 

ambient  conditions  developed  by  Keenan  et  al,8  based  on  the  firing  of  high 
explosives  in  chambers  with  known  vent  areas  and  volumes,  is  given  as 

tg  =  6.28  (Av/V)'-86  (6) 

2/3 

The  above  equation  is  valid  for  A  /V  <0.21.  In  the  current  design  the 
2/3  v 

ratio,  Ay/V  ,  equals  .018,  and  the  duration  time  for  the  quasi-steady 
pressure  is  approximately  1600  ms. 

The  computation  involves  determination  of  peak  residual  overpressure  from 
Equation  (4)  which  when  combined  with  Equations  (5)  and  (6)  yields  the  quasi¬ 
steady  part  of  the  loading  history.  The  junction  between  the  reflected  over¬ 
pressure  history  and  the  quasi-steady  loading  was  smoothed  by  a  curve  inter¬ 
polation  scheme  in  order  to  avoid  a  sharp  transition.  The  resulting  load 
profile  is  shown  in  Figure  3.  This  loading  is  applied  uniformly  at  each 
mesh  point  on  the  inside  wall  in  the  radial  direction  in  the  finite-difference 
structural  response  model  in  the  PETROS  3.5  computer  program9 * * 12  developed  for 
the  BRL.  In  Figure  3  the  load-time  history  inside  the  hemispherical  enclosure 
was  zeroed  out  after  180  ms  to  facilitate  damping  of  small  elastic  oscillations 
and  to  observe  any  residual  deformation  of  the  hemispheric  wall.  The  peak 
reflected  overpressure  was  found  to  be  257.3  kPa,  while  the  peak  residual  over¬ 
pressure  was  approximately  100  kPa. 


IV.  STATIC  STRESS  ANALYSIS 

Prior  to  a  detailed  dynamic  response  study,  a  static  stress  analysis  in 
the  linear-elastic-small  deflection  regime  was  conducted  to  obtain  an  initial 
estimate  of  the  enclosure  wall  thickness.  In  this  investigation  discontinuity 
stresses  at  the  base  of  the  shell  structure  were  ignored.  Since  the  duration 
of  the  reflected  pressure  is  less  than  1.5  ms  compared  to  1600  ms  for  the 
duration  of  the  quasi-steady  overpressure,  an  approximate  static  analysis 
based  on  a  minimum  factor  of  safety  of  2.0  is  considered  to  be  satisfactory. 
For  the  preliminary  investigation,  stress-concentration  near  holes,  cutouts 
and  wall  openings  was  neglected.  However,  the  effect  of  ground-plane  reflection 


8 

W.  A.  Keenan  and  J.  A.  Tamareto,  "Blast  Environment  from  Fully  and  Partially 
Vented  Explosions  in  Cubicles ,"  U.S.  Naval  Civil  Engineering  Laboratory 
Technical  Report  No.  51-027 ,  February  1974. 

9 

S.  D.  Pirotin ,  B.  A.  Berg  and  E.  A.  Witmer,  ,rPETR0S  3.5:  New  Developments 
and  Program  Manual  for  the  Finite-Difference  Calculation  of  Large  Elastic- 
Plastic  Transient  Deformations  of  Multi-Layer  Variable-Thickness  Shells," 
U.S.  Army  Ballistic  Research  Laboratories  Contract  Report  No.  211, 

February  1975  (ADA  007215). 
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Figure  3.  Computed  Pressure-Time  History  Due  to  Internal  Explosive  Blast  Loading  of  the  Hemispherical 
Enclosure. 


of  the  shock  wave  was  included  through  a  load  multiplication  factor  of 
k  =  2.0,  which  in  effect  doubled  the  applied  load. 


To  contain  the  initial  pressure  pulse  in  an  elastic  manner,  only  the  peak 
reflected  overpressure,  P  ,  was  included  in  the  calculation  of  stresses  and 

deflections.  An  equivalent  static  meridional  stress,  o,  can  be  calculated 
from  Figure  2  by  equating  the  resultant  upward  force  due  to  internal 
pressurization  to  the  net  downward  restraining  force  due  to  the  stress 
developed  at  the  clamped  edge,  resulting  in 


a 


RkP 

m 

2h 


(7) 


where 


R  =  8.987  m,  the  median  radius 
k  =  load  multiplication  factor 

However, for  an  assumed  factor  of  safety  of  2.0, a  =  h  a^,  where  is  the 

static  yield  stress.  Substituting  this  value  of  a  in  Equation  (7)  and 
rearranging  terms  results  in  an  expression  for  the  estimated  thickness,  h,  in 
the  form 


h  = 


RkP 

m 


a 


y 


(8) 


The  yield  stress,  a  ,  for  the  wall  material  which  is  1020  steel  is  241.3  MPa. 
Hence  the  wall  thickness,  h,  from  Equation  (8)  is  found  to  be  19  mm. 


Up  to  this  point  no  consideration  has  been  given  to  the  possibility  that 
fragment -induced  damage  to  a  shell  might  result  in  catastrophic  rupture  when 
the  blast  loading  is  applied.  One  should  estimate  the  material  removal  pro¬ 
duced  by  the  impact  of  the  worst  threat  fragment  and  perform  a  local  three- 
dimensional  analysis  of  the  stress  field  to  determine  whether  a  crack  would  be 
propagated  under  such  loading.  This  problem  in  fracture  mechanics  is  difficult 
to  analyze  and  can  be  at  least  partially  circumvented  by  a  conservative 
selection  of  wall  thickness  under  the  assumption  that  the  residual  thickness 
is  capable  of  withstanding  the  peak  quasi-steady  pressure  even  when  a  50% 
depth  of  penetration  has  been  achieved  by  a  part-through  fragment.  The  final 
thickness  chosen  was  25.4  mm  (1  in),  a  material  thickness  which  is  readily 
available.  The  25.4  mm  thickness  corresponds  to  a  stress  level  of  45.5  MPa 
which  when  compared  to  the  yield  stress  results  in  a  final  margin  of  safety 
of  4.3  which  is  satisfactory. 


The  peak  radial  deflection  AR  at  the  pole  is  estimated  from  Reference  10 
as 


If)  •  •  «  • 

R.  <7.  Roark,  "Formulas  for  Stress  and  Strain ,”  Fifth  Edition,  McGraw-Hill 

Book  Company,  Hew  York,  NY,  pp.  96,  451,  1975. 
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AR  * 


(9) 
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R2kPr(l-v) 

2Eh 

where  E,  v  are  elastic  material  properties. 

The  peak  radial  elastic  deflection  at  the  pole  from  Equation  (9)  was 
1.1  mm,  which  is  quite  small  and  is  consistent  with  the  design  objective.  The 
gross  mass  of  the  hemispherical  enclosure  was  approximately  96400  kg  based  on 
a  25.4  mm  wall  thickness.  In  this  study  allowance  was  made  for  the  weight  of 
flanged  material  at  the  base  but  not  for  extra  weight  associated  with  access 
provisions,  welds  or  foundations. 

To  detect  resonance  due  to  coupling  of  the  duration  time  of  the  pressure 
pulse  with  the  natural  vibration  period,  the  time  period,  T,  was  calculated 
from  Reference  10  as 


T  = 


ttR 


V2p(l-v) 


(10) 


where  p  is  the  mass  density.  Further  check  of  interaction  of  the  reflected 
pressure  pulse  due  to  ground  plane  reflection  with  the  elastic  oscillation 
of  the  pole  did  not  reveal  any  significant  problem. 


V.  OPTIMIZATION  STUDY 

An  optimization  study  based  on  equivalent  strength  showed  substantial 
weight  saving  for  a  hemispherical  configuration  at  or  below  6  m  radius  but 
marginal  savings  at  higher  radius  up  to  9  m  due  to  compensating  thickness 
increases.  An  equation  proposed  by  R.  Karpp  et  al  11  for  the  minimum  amount 
of  vessel  material  to  contain  a  specified  charge  is  given  as 


V  = 
m 


(11) 


where  ey  is  the  yield  point  strain  of  the  vessel  material  in  biaxial  tension, 

W  is  the  charge  weight,  p  is  the  density  of  the  vessel  material  and  K  is  an 

-6  a 

empirical  curve-fit  constant  found  to  be  4.08  x  10  m3/kg.  Although  the 
minimum  amount  of  vessel  material  to  contain  a  specific  charge  is  not  the 
governing  design  criterion,  there  may  be  some  interest  in  determining  the 
optimized  value.  If  the  volume  of  vessel  material  is  plotted  as  a  function  of 
the  radius -to-thickness  ratio  of  the  container  as  given  in  Equation  (11),  a 
slow  variation  is  observed  in  the  amount  of  vessel  material  required  to  contain 
the  dynamic  load  as  a  function  of  the  radius-to-thickness  ratio  R/h.  The 
variation  in  material  volume  over  the  design  range  of  350  <  R/h  <  200  is  only 


R.  Karpp ,  T.  A.  Duffey  and  T.  R.  Neal,  " Response  of  Containment  Vessels  to 
Explosive  Blast  Loading ,"  Report  No.  LA-8082,  UC  38,  Los  Alamos  Scientific 
Laboratory ,  Los  Alamos,  New  Mexico,  June  1980. 
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about  3%.  Thus,  the  amount  of  material  required  to  contain  a  specified 
charge  in  this  range  of  configurations  is  essentially  constant.  Very  thin 
wall,  large-radius  vessels  with  R/h  >  400  make  inefficient  use  of  material, 
at  least  for  blast  wave  containment.  On  the  other  hand,  for  thick  wall,  small- 
radius  vessels  with  R/h  <  150,  at  least  12%  or  higher  saving  in  material  vol¬ 
ume  can  be  realized  with  judicious  choice  of  reinforcement  in  critical  sections, 
Unfortunately,  substantial  saving  in  material  could  not  be  achieved  due  to 
constraints  of  minimum  work  space  and  equipment  access  considerations  and  the 
additional  requirement  of  part-through  fragment  containment  with  50%  depth  of 
penetration. 

The  analysis  so  far  applies  only  to  the  containment  of  the  initial  pres¬ 
sure  pulse.  However,  for  long-term  containment  the  volume  of  vessel  material  V 
required  to  contain  the  static  pressure  elastically  can  be  estimated  from  the 
semiempirical  relationship  given  in  Reference  11  as 


v  =  C  ^ 
s  P 

o 


where  C  is  a  constant  with  a  value  of  about  1.3  m3  MPa/kg  for  most  solid 
explosives  and  Pq  is  the  peak  static  pressure.  The  material  volume  appears 

to  be  independent  of  the  radius-to-thickness  ratio  if  the  internal  radius  is 
approximated  by  the  average  radius  of  the  vessel  and  the  usual  formula  for 
equilibrium  of  a  thin  shell  is  used.  Based  on  Equation  (12)  the  material  vol¬ 
ume  required  to  contain  the  static  load  was  found  to  be  approximately  20%  of 
that  required  to  contain  the  initial  dynamic  load. 

A  comparison  of  the  9  m  hemispherical  structure  with  an  equivalent 
9  m  *  9  m  *  7.3  m  rectangular  parallelepiped  all  welded  depleted  uranium  (DU) 
containment  structure12  with  a  .07  m  thick  steel  armor  wall  liner  and  a 
.064  m  thick  roof  liner  indicated  an  increase  of  64  times  in  containment 
capacity  of  equivalent  charge  weight  for  the  hemispherical  structure  with  a 
50%  reduction  in  weight  and  concurrent  doubling  of  the  internal  volume  capacity 
without  any  significant  sacrifice  in  the  minimum  margin  of  safety.  In  addition 
the  simplistic  design  of  the  hemispherical  enclosure,  although  somewhat  diffi¬ 
cult  to  fabricate,  was  a  significant  improvement  for  static  and  dynamic  load 
bearing  considerations.  The  down  time  for  duration  of  residual  overpressure 
was  decreased  substantially  due  to  availability  of  larger  entrance  hole  diam¬ 
eter  and  vent  area  from  3.33  sec  total  venting  time  at  the  R-14  range  to 
approximately  1.5  sec  at  the  R-9  firing  range. 


DYNAMIC  RESPONSE  ANALYSIS 


Response  of  the  structure  subjected  to  transient  loads  from  an  internal 
blast  shown  in  Figure  2  was  conducted  using  the  BRL  version  of  the  PETROS  3.5 
computer  program  [9],  which  employs  the  finite-difference  method  to  solve  the 

_ 

A.  D.  Gupta ,  "Stress  Analysis  of  the  Target  Enclosure  of  the  R-14  Firing 
Range ,  "  Unpublished  Report ,  BRL ,  1980. 
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nonlinear  equations  governing  finite-amplitude  elastoplastic  response  of  thin 
Kirchhoff  shells.  The  model  is  valid  for  large  deflections  and  can  be  em¬ 
ployed  to  treat  the  entire  structure  rather  than  a  small  section. 

A.  Material  Model 


The  uniaxial  tensile  quasi-static  stress-strain  property  for  1020  steel 
which  is  used  as  the  primary  vessel  material  is  shown  by  the  continuous  line 
in  Figure  4.  The  material  is  modeled  in  the  code  as  a  combination  of  three 
linear  segments  indicated  as  the  dashed  curve  in  Figure  4  followed  by  a  per¬ 
fectly  plastic  behavior  and  linear  elastic-plastic  unloading  resulting  in  a 
polygonal  approximation  of  the  experimental  data.  The  strain-hardening  part 
of  the  stress-strain  curve  is  generated  by  a  sublayer  hardening  model  from  a 
weighted  combination  of  elastic  perfectly-plastic  curves,  yielding  a  piecewise 
multilinear  hardening  representation.  Strain-rate  effects  were  neglected, 
which  is  conservative  since  these  effects  increase  the  structural  resistance 
and  thus  reduce  the  total  deformation. 

B.  Finite-Difference  Model 


Since  both  the  responding  structure  and  the  applied  loads  are  symmetric 
with  respect  to  the  vertical  axis  as  shown  earlier  in  Figure  2,  it  is  suffi¬ 
cient  to  model  the  response  of  a  single  pie-shaped  segment  of  the  hemispherical 
enclosure  and  generate  the  entire  structure  by  360°  rotation  of  the  structure 
about  the  axis  of  symmetry,  resulting  in  quite  economical  computer  runs.* 

A  total  of  18  meshes  along  the  surface  and  a  single  layer  through  the 
thickness  were  used  to  represent  the  pie- shaped  segment.  Four  Gaussian  inte¬ 
gration  points  through  the  thickness  were  used  at  each  mesh  for  computational 
purpose.  Total  number  of  mesh  points  did  not  exceed  37.  Initial  configuration 
of  the  finite-difference  model  employed  for  all  subsequent  calculations  is 
shown  in  Figure  5. 


VII.  RESULTS  AND  DISCUSSION 

The  deformed  cross  section  of  the  hemispherical  segment  at  36  ms  relative 
to  the  initial  undeformed  configuration  is  shown  in  Figure  6.  At  this  time 
the  maximum  deflection  occurs  at  the  pole.  The  deflections  are  exaggerated 
due  to  a  high  magnification  factor  of  1000  and  are,  in  face,  small  enough  to 
be  in  the  linear  elastic  range,  in  accordance  with  the  design  objective. 

Figure  7  describes  the  transient  rectangular  components  of  displacement 
in  a  meridional  plane  at  point  A  at  45°  from  the  vertical  axis  of  the  hemi¬ 
sphere.  The  maximum  displacement  at  this  point  is  only  1.028  mm,  essentially 
radially  outward.  Displacements  at  other  locations  are  correspondingly  small 
except  in  the  neighborhood  of  the  pole  of  the  hemisphere  where  a  peak  deflection 
of  1.17  mm  is  observed  at  approximately  36  ms  as  illustrated  in  Figure  8. 
However,  this  displacement  is  less  than  4%  of  the  shell  thickness  so  that 


Since  design  parameters  for  the  entrance  holet  exhaust  openings  and  the 
equipment  access  door  were  not  finalized  in  the  preliminary  concept ,  a 
continuous  configuration  was  considered  for  this  investigation. 
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Figure  4.  Stress-Strain  Property  Modeling 


Figure  6.  Deformed  Configuration  at  36  ms  Corresponding  to  Cycle  No.  1500. 
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TIME  (ms) 

Figure  7.  Transient  Displacement  Component 


geometric  nonlinearities  are  insignificant.  The  larger  response  at  the  pole 
is  attributed  to  focusing  of  flexural  vibratory  energy.2 

The  PETROS  3.5  code  was  run  for  8000  cycles  (192  ms)  in  an  undamped  mode, 
after  which  artificial  damping  was  introduced  to  suppress  the  elastic  oscilla¬ 
tions  which  were  positively  biased  due  to  residual  internal  pressure.  Damping 
was  facilitated  by  zeroing  out  the  internal  pressure.  The  fully  damped  con¬ 
dition  was  achieved  at  cycle  8235  (198  ms)  when  the  final  configuration  was 
found  to  be  identical  to  the  undeformed  configuration  in  Figure  5  with  no 
evidence  of  permanent  plastic  deformation. 

Energy  balance  studies  using  the  code  confirmed  absence  of  plastic  work 
and  numerical  instability.  Both  total  and  kinetic  energies  were  bounded.  The 
fluctuations  of  kinetic  energy  appeared  to  have  twice  the  frequency  of  the 
work  performed  by  the  internal  blast  pressure. 

Transient  strain  components  on  the  outer  and  inner  surfaces  of  the  hemi¬ 
sphere  at  a  point  near  the  edge  are  shown  in  Figures  9a  and  9b,  respectively. 

The  meridional  strain  components  on  the  inner  and  outer  surfaces  are  almost 
in  phase  initially  but  become  out  of  phase  and  unequal  in  magnitude  with 
increasing  time  signaling  the  buildup  of  some  flexural  deformation.  The 
hemisphere  moves  outward  and  inward,  except  at  the  fixed  boundary,  in  a 
spherically  symmetric  breathing  mode  resulting  in  membrane  strains  upon  which 
the  bending  strains  are  subsequently  superposed  due  to  propagation  of  flexural 
waves  from  the  fixed  boundary  towards  the  pole.  Significant  difference  in 
strains  between  the  outer  and  the  inner  walls  at  the  clamped  edge  could  be 
primarily  attributed  to  domination  of  the  response  by  the  bending  strains.  The 
circumferential  strains  indicated  by  continuous  lines  in  Figures  9a  and  9b  are 
zero  as  expected.  Calculations  for  maximum  meridional  stress  based  on  peak 
strain  results  in  a  stress  level  of  48.26  MPa, which  is  equivalent  to  a  safety 
margin  of  4.0.  As  expected  from  elastic  theory,  peak  strains  occurred  at  the 
fixed  edge,  while  maximum  deflection  occurred  at  the  pole. 

The  variation  of  strain  at  the  inner  wall  with  time  at  a  point  near  the 
pole  is  shown  in  Figure  10.  The  continuous  line  depicts  the  circumferential 
strain,  which  is  in  phase  and  very  similar  to  the  meridional  strain  shown  by 
the  intermittent  line.  The  strains  at  the  outer  wall  near  the  pole  exhibit 
elastic  oscillations  of  approximately  the  same  magnitude  and  duration  as  in 
Figure  10.  This  behavior  indicates  substantial  weakening  of  the  flexural 
wave  near  the  pole  and  domination  of  meridional  and  circumferential  strains 
by  the  membrane  component  of  strain  due  to  elastic  vibration  of  the  wall  in 
the  breathing  mode.  The  peak  meridional  stress  at  the  pole  was  calculated 
based  on  elastic  equations  and  was  found  to  be  approximately  25  MPa,  which  is 
considerably  lower  than  the  maximum  stress  at  the  clamped  edge.  The  stress 
level  is  equivalent  to  a  safety  margin  of  8.6  based  on  the  yield  stress. 

Both  strain  components  are  relieved  completely  upon  damping  at  approxi¬ 
mately  198  ms.  An  isometric  view  of  the  fully  damped  configuration  generated 
by  360°  rotation  of  the  pie-shaped  segment  about  the  axis  of  symmetry  is 
shown  in  Figure  11.  The  view  through  Section  A-A  in  this  figure  depicts  the 
final  configuration  upon  damping  superposed  on  the  initial  configuration.  The 
coincidence  of  the  two  configurations  at  a  high  magnification  ratio  of  1000 
indicates  the  absence  of  any  plastic  deformation  and  confirms  small  strains 
and  deformations  throughout  the  structure  in  accordance  with  earlier  results. 
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Figure  9 
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DAMPING  start 


Surface  Strains  at  the  Inner  Wall  Near  the  Pole 


VIII.  CONCLUDING  REMARKS 


It  has  been  demonstrated,  through  use  of  a  rigorous  nonlinear  shell 
response  methodology,  that  it  is  possible  to  design  a  containment  structure 
with  a  hemispherical  configuration  in  an  efficient  and  cost  effective  manner. 
The  methodology  could  be  easily  extended  to  structural  optimization  studies, 
resulting  in  considerable  cost  savings  provided  internal  volume  and  access 
considerations  could  be  met. 

In  spite  of  simplifying  assumptions  and  limitations  of  the  PETROS  3.5 
version  of  the  shell  response  analysis  code  which  neglects  transverse  shear 
deformation  and  rotatory  inertia,  the  analysis  gives  a  clear  insight  into  the 
initial  loading  mechanism  due  to  structural  resistance  and  interaction  of 
various  components  of  the  response.  However,  an  examination  of  the  character¬ 
istics  of  the  hemispherical  structure  reveals  the  following: 

1.  The  9  m  radius,  .0254  m  thick  hemispherical  enclosure  is  an 
efficient  protective  structure  capable  of  withstanding  internal  blast  from  a 
29.03  kg  TNT  charge  with  assured  structural  integrity. 

2.  The  structure  is  capable  of  successful  containment  of  combustion 
products  and  fragments  with  sufficient  mass  and  velocity  to  achieve  a  50% 
depth  of  penetration  with  a  satisfactory  margin  of  safety. 

3.  Peak  deflection  occurs  at  the  pole  due  to  elastic  oscillations  of  the 
structure  in  the  breathing  mode  resulting  from  focusing  of  vibratory  energy 

at  the  pole. 

4.  Peak  strain  occurs  at  the  clamped  edge  and  exhibits  considerable 
difference  in  strain  magnitude  between  the  inner  and  outer  surfaces  due  to 
bending  waves  originating  in  this  region. 

5.  The  ratio  of  the  vent  area  to  the  internal  volume  is  small  enough  to 
result  in  a  slow  rate  of  venting  and  an  extended  venting  time  of  1600  ms  for 
the  quasi-steady  residual  overpressure  to  blow  down  to  the  external  ambient 
conditions . 

6.  Cumulative  damage  effect  due  to  repeated  test  firings  could  conceiv¬ 
ably  cause  low  cycle  fatigue  of  the  structure,  and  a  periodic  inspection  of 
the  internal  surface  and  joints  for  cracks  in  critical  regions  is  recommended. 

7.  Future  work  should  be  directed  to  modelling  of  the  enclosure  structure 
with  wall  openings  for  the  equipment  and  personnel  access  doors,  protective 
walls  for  X-ray  equipments,  detailed  analysis  of  critical  joints  and  stress 
concentration  due  to  holes  and  cutouts  in  comer  regions. 
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LIST  OF  SYMBOLS 


exponential  decay  parameter 
yield  strain  in  biaxial  tension 

Poisson's  ratio 
density  of  material,  kg/m3 
meridional  stress,  kPa 
static  yield  stress,  kPa 

available  vent  area,  m2 

modulus  of  elasticity,  MPa 

empirical  curve-fit  constant  =  4.08  *  10_e  m3/kg 
peak  reflected  overpressure,  MPa 

peak  quasi -static  pressure,  MPa 

reflected  overpressure,  MPa 

mean  radius  of  the  hemisphere,  m 
peak  radial  deflection  at  the  pole,  m 
time  period  of  oscillation,  ms 
internal  enclosure  volume,  m3 
minimum  volume  of  vessel  material,  m3 

volume  of  vessel  material  to  contain  static  pressure,  m3 
equivalent  charge  weight  of  explosive,  kg 

scaled  distance  of  the  wall  from  the  charge  location,  m/kg1^ 
wall  thickness,  m 
blow-down  time,  ms 

positive  phase  duration  of  impulse,  ms 

elapsed  time  from  arrival  of  blast  wave  at  inner  surface  of  hemisphere, 
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Defense  Advanced  Research 
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ATTN :  Tech  Lib 
1400  Wilson  Boulevard 
Arlington,  VA  22209 

2  Director 

Federal  Emergency  Management 
Agency 

ATTN:  Mr.  George  Sisson/RF-SR 

Technical  Library 
Washington,  DC  20301 
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Defense  Intelligence  Agency 
ATTN:  DT-2/Wpns  &  Sys  Div 
Washington,  DC  20301 

1  Director 

National  Security  Agency 
ATTN:  E.  F.  Butala,  R15 
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Defense  Nuclear  Agency 
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2  Commander 
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ATTN:  FCPR 
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ATTN :  FCPRL 
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US  Army  Air  Mobility  Research 
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Eugene  de  Luca 
Watertown,  MA  02172 


32 


DISTRIBUTION  LIST 


No .  of 

Copies  Organization 


1  Program  Manager 

US  Army  BM>  Program  Office 
ATTN :  John  Shea 
5001  Eisenhower  Avenue 
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US  Army  Engineer  School 
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ATTN:  DRDME-EM,  D.  Frink 

Fort  Belvoir,  VA  22060 
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Dover,  NJ  07801 
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US  Army  Armament  Materiel 
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US  Army  ARRADC0M 
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ATTN :  DRDAR-LCB-TL 
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and  Development  Command 
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ATTN:  ATCA-C0,  Mr.  L.  C.  Pleger 
Fort  Leavenworth,  KS  66027 

1  Commandant 

Interservice  Nuclear  Weapons 
School 

ATTN:  Technical  Library 
Kirtland  AFB ,  NM  87115 
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P.O.  Box  10127 
Denver,  CO  80210 

1  IIT  Research  Institute 
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1  Massachusetts  Institute  of 

Technology 
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ATTN:  Dr.  W.  E.  Baker 
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1  Washington  State  University 
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ATTN:  G.  R.  Fowles 
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1  Black  &  Veatch  Engineers 
and  Architects  Inc. 

ATTN:  W.  V.  Hill 
7611  State  Line  Drive 
Kansas  City,  MO  64114 


Aberdeen  Proving  Ground 

Dir,  USAMSAA 

ATTN :  DRXSY-D 

DRXSY-MP ,  H.  Cohen 
Cdr .  USATECOM 

ATTN:  DRSTE-TO-F 
Dir,  USACSL 

Bldg.  E3516,  EA 

ATTN:  DRDAR-CLB-PA 
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USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports. 

1 .  BRL  Report  Number _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) _ 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name: 


Telephone  Number: 
Organization  Address: 
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- -  -  -  -  - - -  FOLD  HERE  - 

Director 

US  Army  Ballistic  Research  Laboratory 

ATTN:  DRDAR-BLA-S 

Aberdeen  Proving  Ground,  MD  21005 
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Director 

US  Army  Ballistic  Research  Laboratory 
ATTN:  DRDAR-BLA-^ 

Aberdeen  Proving  >  round,  MD  21005 
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